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I. INTRODUCTION 

Coherent states are at the heart of semiclassical de¬ 
scriptions of generic quantum systems and have proven 
to be a versatile tool in a multitude of physical problems. 
In the general literature^— , mainly two types of coher¬ 
ent states are typically distinguished: The first type, the 
coherent states of the harmonic oscillator, was already 
investigated by Schrodinger— shortly after the birth of 
quantum mechanics, while SU(2) coherent states in the 
Hilbert space of a spin of general length S were added in 
the early 19709^. 

Both types of coherent states share a list of wellknown 
properties which constitute the basis for their prominent 
role in semiclassics: (i) The coherent states can be gen¬ 
erated by a unitary transformation from an appropriate 
reference state. In the oscillatory case this state is the 
ground state of an harmonic, while for spins one uses 
the heighest-weight state in some arbitrary basis. As a 
result, the coherent states are (ii) (over-)complete, (iii) 
eigenstates of simple operators generic to the system, and 
(iv) they have minimum uncertainty products with re¬ 
spect to an obvious choice of variables. Moreover, (v) 
coherent states show a coherent time evolution perfectly 
mimicking the classical limit under appropriate Hamilto¬ 
nians. For oscillator coherent states such a Hamiltonian 
is the one of the harmonic oscillator itself, and for the 
spin case the Zeeman Hamiltonian (coupling the spin to 
an external magnetic field) plays an analogous role. 

In the present work we argue that one can extend the 
above list by general statements about correlations and 
fluctuations within coherent states. Specifically we con¬ 
sider the coherent expectation value of a product of two 
arbitrary operators. For such expectation values we de¬ 
rive systematic expansions in the quantum parameters 
h or 1 /S which involve only coherent expectation values 
of single operators and their commutators with the sys¬ 
tem variables. These expansions are a versatile tools for 
the study of the semiclassical regime of generic quantum 
systems. As an important finding, the energy fluctua¬ 
tions of an arbitrary Hamiltonian are generally proven 
to be in leading order entirely due to the time depen¬ 


dence of the classical variables. These results add to the 
above list of properties of coherent states. Reflecting the 
widespread use of the latter objects, we apply our find¬ 
ings to the Lipkin-Meshkov-Glick model originating from 
nuclear physics^, to the Dicke model describing superra¬ 
diance in quantum optics^, and to coherent intertwiners 
of spin networks occurring in the loop approach to quan¬ 
tum gravityi^^. 

This paper is organized as follows: In section [TT] we 
review and summarize important properties of oscilla¬ 
tor coherent states and SU(2) coherent states. The an¬ 
nounced results on the coherent expectation values of ar¬ 
bitrary operator products are derived in section m and 
discussed there on a general footing. Some technical de¬ 
tails of the calculations are deferred to appendix [A] Sec¬ 
tion [IV] contains the application of our general findings to 
the Lipkin-Meshkov-Glick model, and the Dicke model is 
treated in sectionlVl In sectionlVDwe turn to the study of 
coherent intertwiners of spin networks investigated in the 
covariant approach advocated by Loop Quantum Grav¬ 
ity. Here we derive semiclassical corrections to expecta¬ 
tion values in terms of universal expansion coefficients 
depending only on the network geometry. We close with 
a summary and an outlook in section IVIII 

II. COHERENT STATES 

We now briefly review, using standard notation, dis¬ 
tinctive properties of oscillator coherent states and SU(2) 
coherent states. 

A. Coherent Oscillator States 

The harmonic oscillator is described by 

'Hh = ^ (p 2 + wV) = fiw (a + a + 0 (1) 

with 
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fulfilling 

[p, q] = - ^ [a, a + ] = 1. (3) 

l 

The system has an equidistant spectrum labelled by n £ 

{ 0 , 1 , 2 ,...}, 

H h \n) = fiio ^n+ \n). (4) 

Coherent states of the harmonic oscillator are eigenstates 
of the lowering operator a with complex eigenvalues a , 

a\a) = a\a). (5) 


for the direction s = (sindcost/?, sindsin ip, cos'd). For 
generic systems expectation values within these states 
provide a natural approach to the classical limit given by 
H —> 0, S —> oo while fiS is kept constant. 

Introducing the usual basis of eigenstates of S z 
( S z \m) = hm\m)) coherent states can be generated from 
| S') by a unitary rotation, 


\#,<P) = U^, V )\S) 

= (TTW^ 

with^ 


(14) 

(15) 


They are generated from the ground state via 
| a) = exp (aa + — a* a) | 0 ) 

,2 N 

n —0 


= exp [ — — |a| 
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The parameter a is naturally decomposed into its real 
and imaginary part as 



( 8 ) 


U(d,ip) 

= exp HOI 

(sin <pS x — cos 

<ps») J 

(16) 


_ £ zS /H e riS z /h e ~zS + /h 


(17) 


= e~ zS+ / h e 

-rjS z /h^zS /h 


(18) 

and 





z(d,tp) = 

d ■ 

tan — e l{p , 

d 

77 (d) = 2 In cos — . 

(19) 


Expanded in the above basis SU(2) coherent states read 


Denoting an expectation value within a coherent state 
© by (•) it holds 


(q) = £ > (P) = 7T ■ (9) 

Coherent states maintain their shape in the time evolu¬ 
tion of the harmonic oscillator, 

e-* Hht \a) =e-*“ t |ae -i “ t >, (10) 

and the time dependence of the expectation values © 
follows exactly the classical motion of the harmonic os¬ 
cillator. This fact justifies the term ‘coherent states’ and 
relies on the equidistance of the spectrum. The latter 
property is shared by a quantum spin of arbitrary length 
in a magnetic field and leads there to a coherent Larmor 
precession, as we will discuss in section Hi Bl 

Moreover, coherent states minimize uncertainty prod¬ 
ucts, 

ApAq = H/2 (11) 

and fulfill an (over-)completeness relation, 

— / d 2 a\a)(a\ = 1. (12) 

B. SU(2) coherent states 

In the Hilbert space of a spin of length S an SU(2) (or 
spin) coherent state |d, </?} is defined by the equation 

s- 5|d, ip) = hS I'd, ip) (13) 
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z m \m) ( 20 ) 


= E 


m=—S 


2S \ 5 / /d 
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■ ( 21 ) 


The analog of the harmonic oscillator for SU(2) coherent 
states is the Zeeman Hamiltonian 

U z = -S-h ( 22 ) 

coupling the spin to a magnetic field h. The spectrum 
consists of 25 + 1 equidistant energy levels, and the cor¬ 
responding time evolution of SU(2) coherent states is a 
coherent Larmor precession, which is most easily seen 
when putting, without loss of generality, the field direc¬ 
tion along the 2 -axis, 

p) = e~^ Sht \d, p + ht). (23) 

The latter finding is completely analogous to the har¬ 
monic oscillator having a semi-infinite equidistant spec¬ 
trum. 

As further standard properties shared with coherent 
oscillator states, SU(2) coherent states have a minimum 
uncertainty product 


A(ei • S ) A(e *2 ■ S) = 


H 2 S 

2 


(24) 
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with ei,e* 2 ,s being an orthonormal system, and their 
(over-) completeness can be expressed as 


1 = 


25+1 

47T 

25+1 

7T 

25+1 

7r 


/»27T n 7T 

/ dp d#sin$|$, (25) 

Jo Jo 


(i+ki 2 ) ; 


l*X*l 




e zS / s |5)(5|e 2S+ / R 


(1- 


2)2(S+1) 


(26) 

(27) 


where |z) = |$, +). Further below it will be useful to 
change reference state 15) in Eq. (1271) to an arbitrary 
SU(2) coherent state by applying the unitary transfor¬ 
mation given in Eqs. 


i = 


25+1 

7r 

25 + 1 

7r 


U 


' e wS ~/ h \S){S\e^ s+ / h Tr+ 

W (l + | u ,|2)2(S+r) 

l2 e w ' s ~/ h \z){z\e^ + / h 


(1- 


w 


2)2(S+1) 


(28) 


with 5 = USU+. 


III. CORRELATIONS 

We now derive general theorems for the expectation 
values of operator products within coherent states. 


A. Oscillatory Systems 


1. General Correlation Functions 

Let A, B be two operators being functions of the two 
canonical operators p, q (or, equivalently, a, a + ). Using 
the completeness relation (1T51) the expectation value of 
AB within a coherent oscillator state can be formulated 
as 


(a\AB\a) = i J d 2 /3e~W 2 {0\U+AU a U+e? a+ \0) 

■(0p a U a U+BU a \0), (29) 

where U a is the unitary operator on the r.h.s. of Eq. m, 
and 

U+e 0a+ 10) =e-^ a l 2+a/3 e (/3 - a)a+ |0) (30) 


such that 
{a\AB\a) 


J d 2 /3e-l^-“l 2 (0|f/+ AU a e^~ a ^ a+ \0) 

■(0\e(P- &)a U+BU a \0) 

J d 2 /3e-^ 2 (O\e- pa+ U+AU a e 0a+ \O) 

■(0p a U+BU a e-P a \0), (31) 


where we have shifted the integration variable and used 
e -/3“|0) = |0). The remaining operator products can be 
expanded into series of iterated commutators according 
to 


oo 1 

e*y e - Y = £-[*,+]„ (32) 

n =0 

with [A,r] 0 = Y and [X,Y] n = [X,[X,Y] n _ i]. Upon 
performing the integration the two infinite series shrink 
to a single one yielding 


(a\AB\a) = £-(0| [~a+, U+AU a \ JO) 

n=0 

■(0\[a,U+BU a ]j0) (33) 

OO 1 

= \iU a a+U+,A] n \a) 

n —0 

•(a| [iU a aU+,B] n \a) 

oo 1 

= [ ia+ ^ A ]J a ) 

n —0 

■(a\[ia,B] n \a). (34) 

In the last step we took into account that U a a + U+ and 
a + differ just by a constant which commutes with any 
operator. Thus, we have arrived at an expression for the 
expectation value of product of two operators within co¬ 
herent states in terms of a sum over products of such 
coherent-state expectation values which involve only one 
of the operators. An alternative form of the above ex¬ 
pansions can be given via Eq. (ED as 


(a\AB\a) = ^T(0|U+A[/ Q |n) (n\U+BU a \0 ), (35) 

n —0 


which of course just expresses the completeness of the 
states In) and provides an alternative way to derive 

Eq. ED. 

Moreover, using the definition m Eq. (1341) can be 
rewritten as 


(a|Ai3|a) 



n —0 








(36) 


Since each commutation of p, q with A or B yields a 
factor of ft all expectation value on the above r.h.s. are of 
the same order in ft. Thus, Eq. (157)1) is indeed a systematic 
expansion in ft of the coherent-state expectation value of 
an arbitrary product of two operators. The zeroth order 
equals the classical result, and for a general correlation 
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function one has the semiclassical expansion 


Cab ■— {a\AB\a) — {a\A\a){a\B\a) 


OO 


- E 


n= 1 

'a 


77 , 12 " 






(37) 


Choosing A = B we obtain a general expression for the 
variance of an hermitian operator A , 


Relations of the type (TTTT|) . fPT2jt were already found in 
RefJ^ on the example of specific Hamiltonians. The re¬ 
sults here are derived for arbitrary systems and are based 
on the very general expansions mm for correlation 
functions and fluctuations. 

The fact that the system will in its time evolution in 
general not strictly remain in a coherent state, i.e. de¬ 
coherence occurs, is reflected by the higher contributions 
to the energy variance. Indeed, for a harmonic oscillator 
m the time evolution is strictly coherent and we have as 
an identity 


(AH) 2 



n =1 





(38) 


where each term in the semiclassical expansion is non¬ 
negative. 


2. Energy Fluctuations 


Considering A = PL as an Hamiltonian, the correspond¬ 
ing energy fluctuation reads in leading order in H 


(AH) 2 


H 

2 

H 

2 




(39) 


(40) 


(A U h f = ^ 


-(W 


(d t n) 2 


(43) 


for all times t > ti and without any higher correction. 

Finally, it is straightforward to extend the above re¬ 
sults for general operator products to the case of TV > 1 
degrees of freedom; details are sketched in appendix [A] 
For the energy variance one finds in leading order in h 


N 


(AH) 2 = -E 


|lv^.,Hl) +(j 


Pa 

\[CT a 


1 \ 2 


,n 


+o ( h 2 ) 


£ N 
ri \ 

9 2 ^ 


w 0 (dtq a ) 


(■ dtPa} 2 
Wo 


(44) 

(45) 


with operator pairs p a , q a and frequencies w a , and the 
analog of Eq. (H21) reads 


where we have replaced, according to the Heisenberg 
equations of motion, the commutators with time deriva¬ 
tives. Indeed, if the system is prepared at some initial 
time t = ti in a coherent state we have (cf. Eq. ()8|)) 


(AH) 2 


h 

2 


N 


E 


Wo (<9t£o) 2 + 


(d t n a ) 2 

Wo 


+ O (n 2 ) . (46) 


(d t q)=d t £ , (d t p)=d t Tr (41) 


B. Spin Systems 


and 

w (c> t £) 2 + + O ( h 2 ) (42) 

at t = ti. In the subsequent time evolution governed by 
the Hamiltonian PL the state of the system will, for not 
too large times, approximately be coherent with time- 
dependent parameters £(f), 7r(t) playing the approximate 
role of classical Hamiltonian variables. Thus, in this 
semiclassical regime the fact that a coherent state has 
a finite energy variance, i.e. it is not an eigenstate of 
the Hamiltonian, is in leading order in h just expressed 
by the fact that the classical Hamiltonian variables have 
a nontrivial time dependence, i.e. the system is moving. 
This result complements the historical Ehrcnfest theorem 
stating that expectation values of observables follow the 
classical equations. 


1. General Correlation Functions 

We consider again two arbitrary operators A, B which 
are now functions of a spin operator S. The expectation 
value of the product AB within an SU(2) coherent state 
| z) for spin length S can be formulated as 

= ^/ u + m 

■{z\e- w§ ~ /h Ae w§ ~ ' h \z) 
■(z^^^Be-^/^z), (47) 

where we have used the completeness relation in the form 
(l28ll and the observation 

e -ws+/h\z) = f/ e -^ s+ / R |0) = | z) (48) 
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with U given in Eqs. m-m- Employing now again the 
expansion (1321) and performing the integration leads to 


(z\AB\z) 


y, (2S-n)! / 

h "W v 




The above equation is the spin analog of the result (1551) . 
Again all iterated commutators are of the same order in 
h and 5 whereas the prefactor of the n-th term carries a 
product 25(25 — 1) • • • ( 2S — n + 1) in its denominator. 
Thus, Eq. (1491) is essentially an expansion in the quantum 
parameter 1/5. Note that the spin components S x , S v 
represent the direction perpendicular to the spin polar¬ 
ization of the coherent state | z). Alternatively, the result 
(09| can be written as 


(z\AB\z) 


2 S 


E 

71=0 


25 


E 

71 = 0 


(25 — n)\ / 
n!(25)! \ 



(S\U + AU\ 5 
■ (5 — n 


l -S-,U + AU 

a 

S+,U+BU 


n 



— n) 

\U+BU\S) . 


5 

(50) 

(51) 


Analogously to Eq. (1551) for oscillatory systems, the last 
formulation is just the completeness relation for the 
states | m) and allows for an alternative derivation of the 
central result (1491) . Using the latter, arbitrary correlation 
functions within SU(2) coherent states can be expressed 
in full analogy to Eq. (1371) . 


2. Fluctuations 

For the variance of an hermitian operator A we have 
2 S M , r • i \ 2 


(^) 2 = E 


71=1 


(25 -n)! 
i!(25)! 


-S~,A 


(52) 


The expectation values occurring in leading order can be 
rewritten as 


iS~,A 


= E 


1=1 

3 


iS l ,A 


(53) 


= EI<*I[«MI*>1 > ( 54 ) 


1=1 

where we have observed that | z) is an eigenstate of S z , 
and that S and S are related by an orthogonal matrix, 

3 

5 - (55) 

j=i 


Thus, we have 


(AA) 2 





, (56) 


and by a slight generalization of the above arguments one 
finds for the expectation value of a product of commuting 
operators A, B 


(z\AB\z) = \{z\AB + BA\z) 


= (z\A\z)(z\B\z) 

3 

z 


—E 

2 


i= 1 


^S\A 

a 


I st > B 


+° ( ^2 I ' 


(57) 


The requirement here for a symmetrized operator prod¬ 
uct stems from the fact that for an identity analogous to 
Eq. (1551) to hold products of expectation values involving 
both S x and S v should drop out. 

Choosing now in Eq. (1551) A to be the Hamiltonian 
FL of the underlying system we can write by the same 
arguments as for Eq. (SO]) 


^ = h( d ‘ s ) 2+ °{h) (58) 


with (•) = (z\ ■ |z). To this result the same comments 
apply as to its oscillatory counterpart Eq. (HU1) : If the 
system is initially in an SU(2) coherent state it holds (cf. 
Eq. (USD) 


(d t S^ = HSdtS (59) 

and 

(AH) 2 = (hSf (E ( 5iS -f + O (E) ) (60) 

at initial time t = tj, and for not too large times t > ti 
the system will approximately remain coherent in its time 
evolution under H with s(t) being a classical vector. Our 
finding (1601) is again a manifestation of our previous re¬ 
sult (l42|) : In leading order in the quantum parameter (h 
or 1/5) the variance of the energy is due to the classical 
motion of the system. Findings of the type (l60l) were 
also obtained previously in Ref.— on the example of spe¬ 
cific Hamiltonians. Here we provide a generalization to 
arbitrary systems based on the very general expansions 
(|49|) . (152l) for correlation functions and fluctuations. 

Decoherence effects, i.e. deviations from the coherent 
state with time-dependent parameters s(t) are again in¬ 
dicated by the higher-order terms in the energy variance, 
as we shall investigate on a specific example in section 
IIV Al Conversely, the Zeeman Hamiltonian (l22l) gener¬ 
ates a strictly coherent time evolution with 


(A H z ?= { ^f(d t s) 2 


25 


(61) 
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as an identity for arbitrary times t > ij. 

Similarly as for oscillator systems, the above results 
for general operator products are easily generalized to 
the situation of N > 1 spins of various lengths; details 
can be found in appendix m The leading order of the 
energy variance is given by 


N 


(Anf = J2 


a— 1 


2S7 (8A 




(62) 


and Eq. ([501) is generalized to 


N 


(An ) 2 = £ 


(n 5 a) 2 ^i (at ' a)2 + °(i 


(63) 


IV. THE LIPKIN-MESHKOV-GLICK MODEL 

The Lipkin-Meshkov-Glick (LMG) model is an approx¬ 
imate description of N interacting spin-1/2 systems and 
was originally inspired by nuclear physics ^ 15 ! 16 . More re¬ 
cently this model has been argued to describe two-mode 
Bose-Einstein condensates^ - —, phase transitions in op¬ 
tical cavity QED21 - — , and molecular magnets— More¬ 
over it has been employed to model a spin bat h—' — and 
in studies of quenched dynamics— In the last decade 
a flurry of publications investigating various aspects of 
the LMG model has appeared; as an entry point to the 
recent literature we refer to Refs— - —. 

Concentrating on the sector of maximal spin S = N/2, 
the LMG Hamiltonian reads 

n = ~ hSZ ~ i {lxSXSX + 7 y svsv) > ( 64 ) 

where h can be interpreted as a magnetic field coupling to 
the z-component of the spin while 7x, 7 y parametrize an 
anisotropic interaction among the perpendicular compo¬ 
nents. The factor hS in the denominator is a convention 
common to the literature and leads to a linear scaling of 
energies as a function of S 1. The expectation value 
with in an SU(2) coherent state is given by (neglecting a 
constant contribution) 

(n) = hs(^— h cos’d 

+-y sin 2 d cos 2 p + ^ sin 2 d sin 2 p'j (65) 

and equals the classical energy expression up the renor¬ 
malized parameters 7,: = 7* (1 —1/(25)). Taking coherent 
expectation values of both sides of the Heisenberg equa¬ 
tions of motion one obtains the (semi-)classical equations 


For the energy variance one finds by a direct (and some¬ 
what tedious) calculation 

(An ) 2 = fti + n 2 (69) 


with 


^ = {hS)2 is 


h 2 sin 2 d — 2h ; y x cos d sin 2 d cos 2 p 


—2hj y cos d sin 2 d sin 2 p 

+7 2 (sin 4 d cos 2 p sin 2 p + cos 2 d sin 2 d sin 2 p) 
+7 2 (sin 4 d cos 2 p sin 2 p + cos 2 d sin 2 d cos 2 p) 

—27 x 7 y sin 4 d cos 2 p sin 2 p 


(70) 


being of leading order 1/5 while the contributions sum¬ 
marized in 


n 2 = (hsy 


8 s 2 


1 - 


2 s 


-47x7 y cos 2 d 


+ (7a: (l — sin 2 d cos 2 p) + 'jy (l — sin 2 d sin 2 </?))' 


(71) 


are of order 1/5 2 and higher. Using now Eqs. (I66ll - (l68l) 
we can identify the leading contribution to the energy 
uncertainty as 


Hi 



(72) 


in accordance with the general result (16U1) . The sublead¬ 
ing contributions indicate decoherence effects, i.e. de¬ 
partures from the submanifold of the coherent states in 
the Hamiltonian time evolution, as we now discuss ex¬ 
plicitly on the example of the isotropic LMG model. 


A. The isotropic case 

Putting 7x = 7y =: 7 the Hamiltonian becomes diago¬ 
nal in the states | m) with eigenvalues 

Em/h = ~hm + ^m 2 - ^(5+ 1). (73) 

This eigensystem is simple enough to analytically com¬ 
pute the exact time evolution of coherent expectation 
values ( S(t )): Due to symmetry, the z-component is con¬ 
stant, 


{S z (t)) = hS cosd (74) 

while for the perpendicular components one finds 


ds x 

dt 

ds y 

dt 

ds z 

dt 


h sin d sin p — J v cos d sin d sin p , (66) 

— h sin d cosp + cos d sin d cosp , (67) 
— (lx — ly) sin 2 $ cos 7) sin 7). (68) 


(S+(t)) = hSsmde i ^- ih -^ 1 - 1/{2S))cosi,)t) 


■ (e z=L ^ ±t cos(^i)+icos0sin(^i) ) (75) 

The above closed result relies on the fact that S + couples 
only eigenstates with neighboring indices such that all 
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occurring energy differences are, apart from a constant 
term, linear in to. The first line in Eq. (USD describes a 
classical rotation of the spin according to Eqs. m-m 
whereas the second line contains quantum effects: The 
“spin length” 


oscillator and an SU(2) coherent state reads 

(H) = Hu\a\ 2 + klHS cos d 

+ . HSsmdcosip (a + a) , (80) 

V 2 S 


|(S ,+ (t))| = HS siwd (l - sin 2 t? sin 2 ^ ( 76 ) 

composed from the perpendicular components breathes 
sinusoidally in time. Quantum (quasi-)revivals occur at 
times at t = 2nkS/~/ for any integer k where the state re¬ 
turns precisely to the submanifold of the coherent states. 
These times are large in the semiclassical regime as they 
are proportional to S. 

Regarding small times, we define t =: \fSr and con¬ 
sider the regime yr <C y/S, such that for large S' 1 it 
follows 


\(s + m 

hS sin d 



(jt sin$) 2 /4\ 

S ) 


e ~(lT sini9) 2 /4 


5—1/2 


(77) 


which perfectly matches the classical expression. The 
(semi-)classical equations of motion can be obtained 
analogously as Eqs. HMD-®, 


da 

dt 

ds x 

dt 

ds v 

dt 

ds z 

dt 


— ^ A , 

zua + — r _ HS smi7cos</>, 


frV2S 

—S2 sin $ sin </;, 


H sin d cos ip - -= cos d (a + a) , 

V 2 S 

\ 

sin d sin ip (a + a) , 


V2S 


(81) 

(82) 

(83) 

(84) 


and a direct (but again quite lengthy) calculation of the 
energy variance yields 

(AH) 2 = + 0 2 (85) 


i.e. the spin expectation value ( S(t )) shows a gaussian 
decay with time scale At = y/2S / (7 sin d). On this time 
scale, sometimes known as Ehrenfest time^S, departures 
between classical and quantum dynamics become sizable. 
The above finding for At is consistent with a heuristic 
uncertainty argument in the following sense: Replacing in 
AH At > H the energy uncertainty with y/Ti^ one obtains 
a lower bound for At being proportional to HS which is 
a constant independent of S in the semiclassical regime. 
Thus, this lower bound is consistent with the above result 
which grows with the square root of S. 


V. THE DICKE MODEL 

The Dicke model describes the superradiant interac¬ 
tion of a single cavity mode of a radiation field with N 
two level systems (atoms)^. Although introduced al¬ 
ready in the 1950s, this model continues to be investi¬ 
gated under various aspects; as a guide to the recent 
literature see e.g. Refs .— —. 

Focusing again on the sector of maximal spin S = N/ 2, 
the Dicke Hamiltonian can be formulated as 

H = Hua+a + nS z + ^=S X (a++ a) (78) 
v2 S 

= - (p 2 + u 2 q 2 ) + flS~ + \^Jj^S x q , (79) 

where the parameters w, H, and A have all dimension 
of inverse time. In the classical limit, the superradiant 
phase, characterized by a finite bosonic occupation in 
the ground state, occurs for A 2 > fiw. The expectation 
value of the Hamiltonian within a tensor product of an 


with the leading-order term 

Hi = (hu>) 2 \a\ 2 + ^77-n 2 sin 2 d 

ZO 


+ < fis 4 


— (sin 2 d sin 2 ip + cos 2 d) (a + a) 2 
Zo 


+ sin 2 cos 2 ip 


+Hu—j^HS sindsin^ (a + a) 

V 2 S 

A (HS) 2 

-V. — ——— cos d sin d cos ip (a + a) 
v 2 S S 

and the subleading contributions 
A 2 


(2 sy 


(<S"> 2 + <S‘> 2 ) . 


Finally, comparison with Eqs. (l8lT) - (l84l) shows 


( 86 ) 


H 2 = (HS) 2 (sin 2 d sin 2 <p + cos 2 d) (87) 

(2S) 


( 88 ) 


^1 — 77 w (<%£) + 


2 . (dt^) 2 \ , (HS) 2 


w 


2 S 


(d t s) 2 , (89) 


in accordance with the general results (H51) . (1551) . and 
(lATl) . 

Note that the higher-order contributions (1551) describ¬ 
ing decoherence effects depend only on the coupling pa¬ 
rameter A but not on the frequencies w, H, in accordance 
with the fact that spin and oscillator show perfectly co¬ 
herent time evolutions in the absence of coupling. An¬ 
other distinctive feature of the result (1551) (compared to 


















e.g. Eq. (ED ) is its simplicity which calls for further ap¬ 
plications. In fact, an extensive numerical study of the 
dynamics of the Dicke model in the semiclassical regime 
was performed recently in Ref42. Here the initial condi¬ 
tion was, as in the present work, a tensor product of an 
an oscillator and a spin coherent state, and it is straight¬ 
forward to evaluate the above H2 in terms of such dynam¬ 
ical data. In particular, it is an interesting speculation 
whether or not Q 2 behaves differently in the regular ver¬ 
sus (quantum) chaotic regime as studied in Ref4£. An¬ 
other aspect is to compare the Ehrenfest times At found 
numerically with estimates according to y/Tl^At > h. 


VI. COHERENT INTERTWINERS IN SPIN 
NETWORKS 


We now apply our general findings on coherent expec¬ 
tation values of operator products to spin network states 
as studied in Loop Quantum Gravity (LQG)^ ’ 44 ’ 45 . In 
brief, a spin network is a collection of points (called ver¬ 
tices or nodes ) in (typically) three-dimensional space con¬ 
nected by one-dimensional curves (edges). Each edge is 
assigned a spin of individual length, and a spin network 
state in the tensor product of all those SU(2) represen¬ 
tations is defined by the additional requirement that all 
spins joining in a given node are coupled to a total singlet. 
The latter property implements the Gauss constraint on 
the holonomy and flux variables used in LQ C 46 ’ 47 . 

A convenient parametrization of spin network states 
are coherent intertwiners as introduced by Livine and 
Spezialoii. Fixing an TV-valent node (connecting N 
edges), one considers a tensor product 

N 

l*>:=0K,¥>a> (90) 

a— 1 


of SU(2) coherent states describing the spin on each edge. 
A coherent intertwiner is then defined by the projection 
of this object onto the singlet subspaceii 


\$)s 




(91) 


where the denominator takes care of the normalization. 
The projection operator can be formalized by a Haar in¬ 
tegration over all uniform rotations of the N spins (group 
averaging), 



(as a factor of h will occur below in the Planck length 
squared). In particular, a coherent intertwiner is by con¬ 
struction invariant under arbitrary rotations of all spins 
meaning 

(e*) I*>. = °- ( 94 ) 

Moreover, nodes in a spin network allow for a geomet¬ 
ric interpretation in terms of convex polyhedral. From 
a classical point of view, this relies on a theorem due to 
Minkowski^. It states that given N unit vectors s a and 
N positive numbers A a fulfilling ^aS a = 9, there is 
a unique convex polyhedron with N faces such that s a 
is the normal to the a-th face and A a is its area. Thus 
choosing as areas the quantum numbers S a , the classical 
closure relation 

N 

s *Sa = 0 (95) 

< 2=1 

ensures that the geometric information contained in the 
state (TO1) encodes a convex polyhedron. The quantum 
counterpart of the relation (liJ5l) is equation (liTH) giving 
rise to the notion of a quantum polyhedron In the 
framework of LQG, the spin operators representing the 
faces of the polyhedron are, up to a prefactor, considered 
to be flux operator s 46 ’ 47 

E a = 87 n £ 2 pS a (96) 


with 7 being the Immirzi parameter and the squared 
Planck length £ 2 P = HG/c 3 . 

Let us now explore expectation values within coherent 
intertwiners. Here one can concentrate without loss of 
generality on operators unchanged by uniform rotations 
since for any operator being the sum of a rotationally 
invariant part and terms without this property, only the 
former will contribute. Any rotationally invariant op¬ 
erator Q commutes with the projector onto the singlet 
space, [Q,P] = 0, such that PQP = QP = PQ. There¬ 
fore we can use the result ED to obtain a semiclassical 
approximation to the expectation value within a coherent 
intertwiner, 


= <*|Q|*> 


($|PQ + QP|$) 

2($|P|$) 


(97) 


+E^E( $ l[ iS «>G] I *)<£(*)+ ■■• 

a=1 i=l 


(98) 


with 






(99) 


with n = (sin d cos y?, sin d sin yj, cos d). Here and in what Thus, the expectation value of Q is in leading order just 
follows we take all spin operators to be dimensionless given by the expectation value of the unprojected state 
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and the normalization factor in the definition m 
drops out. For the subleading corrections one needs to 
determine the coefficients (1551) . Here both numerator 
and denominator are conveniently formulated in terms 
of Haar integrations as shown explicitly in Eq. (1531) . In 
the semiclassical regime studied here where all spins are 
long, VaS'a 3> 1, these integrals become amenable to a 
saddle-point approximation as worked out in Refi 11 . For 
the denominator one finds for a general IV-valent node 


<$|P| $) 


1 

s/ 7r det if 


tr (if- 1 ) 
4V 7T det if 


( 100 ) 


where 


N 

H^ = J2 S a{S ij -si4) ( 101 ) 

a= 1 


is twice the negative Hessian of the saddle point expres¬ 
sion, and the details of the calculation can be found in 
appendix [Bj Since if is a linear combination of geomet¬ 
ric projection operators (<5 lJ — s l a s J a ) with positive coeffi¬ 
cients, its eigenvalues are nonnegative, and zero eigenval¬ 
ues only occur in the degenerate case where all vectors 
s a are collinear, which we shall not consider here. Thus, 
the eigenvalues of if can be taken to be positive, and the 
determinant can be formulated more explicitly as^i 

det if = ^ y S a S b (s a x s b ) 2 

ab 

-l^S a S b S c \(s a xs b )-s c \ 2 (102) 

abc 

with T = ^2 a S a ■ The semiclassical limit of a quantum 
polyhedron is obtained by rescaling all quantum numbers 
as S a i->- A S a with some integer A> 1. Thus the leading 
term in Eq. (HOOD (already obtained in Ref—) is of order 
A -3 / 2 while the subleading correction scales like A -5 / 2 . 
The numerator in Eq. (1551) can be evaluated via saddle 
point approximation in a similar fashion (see appendix 
HD giving, again for a general TV-valent node fulfilling the 
classical closure relation (1551) . 


$ 


iS a ,P 


$)=S a 


S a X (if 1 S 0 ) 
VirdetH 


(103) 


such that for the coefficients themselves we have the 
amazingly simple result 


us that the coherent-intertwiner expectation value of any 
(rotationally invariant) operator is in leading order given 
by the expectation value of the unprojected tensor prod¬ 
uct of SU(2) coherent states, and the leading correction 
scales with the inverse of the spin lengths. 

The coefficients (11041) are universal in the sense that 
they are the same for any operator Q. Making use of the 
symmetry of if they can also be formulated as 

cim = s a j 2 z ijk (h'T ( 105 ) 

jH 

implying the sum rule 

N 

^Cl(<F) = ~J2 eijk {H~ 1 ) kl H lj = 0 , (106) 

a— 1 jkl 

which also follows from the definition (1551) and the quan¬ 
tum closure relation (1551) . Thus, the sum rule (11061) holds 
independently of the fulfillment of the classical closure 
relation (l95l) which underlies the explicit result (11041) . 

Moreover, it is interesting to note that the matrix if 
can be interpreted as the inertia tensor of a distribution 
of masses S a whose positions are given by the unit vec¬ 
tors s a . By the same token, the classical closure relation 
®>D states that the center of mass of this distribution lies 
in the origin of the chosen coordinate system. In partic¬ 
ular, if is proportional to the unit matrix (such that the 
expansion coefficients (1551) vanish) if the node has the 
shape of an archimedian body such as a regular tetrahe¬ 
dron. We leave it to further studies to explore further 
possible consequences of the above analogy. 

Very typical examples of rotationally invariant opera¬ 
tors are volume operators of polyhedr a 48 ' 50 ! 51 . The sim¬ 
plest nontrivial case of a quantum polyhedron is given 
by a tetrahedron, i.e. a 4-valent node— The volume 
operator can be formulated as 

V =^yf\E 1 -(E 2 x E 3 )\ (107) 

using any three of the four flux operators. Squaring this 
expressions and striping all prefactors one is led to con¬ 
sider the expression 

Q = S 1 -(S 2 x S 3 ) . (108) 

acting on the Hilbert space defined by the constraint (1941) . 
The study of this operator in the semiclassical limit has 
attracted quite a deal of interest recentl y 48 ’ 53- — . For the 
expectation value within coherent intertwiners one finds 
from Eqs. (EED-CESS 


d($) = SaSa X (H-'sa) . (104) 

The expression (11031) is of order A -3 / 2 while the coeffi¬ 
cients (11041) are independent of A and vanish if the matrix 
H is proportional to the unit matrix. Thus, polyhedra 
where all eigenvalues of if are degenerate enjoy an en¬ 
hanced classical character in the sense that the leading 
order of semiclassical corrections to general expectation 
values (HJ51) vanishes. In the general case, Eq. (I551l tells 


,<$|Si(S 2 x S 3 )|*> s = SiS 2 S 3 


Si (s 2 X S3) 


+ “ (s*l X (s 2 X S 3 ) • (si X (if 4 Si)) + C.p.) 

+ ••• (109) 

As before, the form of the subleading corrections here 
holds if the classical closure relation (1551) is fulfilled. 
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VII. SUMMARY AND OUTLOOK 


Appendix A: N > 1 degrees of freedom 


We have derived general systematic expansions with 
respect to quantum parameters of expectation values of 
products of arbitrary operators within both oscillator co¬ 
herent states and SU(2) coherent states. These results 
are a versatile tools for the study of the semiclassical 
regime of generic quantum systems. In particular, we 
prove that the energy fluctuations of an arbitrary Hamil¬ 
tonian are in leading order entirely due to the time de¬ 
pendence of the classical variables, a result very general 
and very intuitive at the same time. 

Our findings offer many possibilities for application in 
various fields of physics. Here we have specifically stud¬ 
ied the Dicke model stemming from quantum optics, and 
the LMG model originating from nuclear physics. For 
the latter system we have investigated decoherence ef¬ 
fects (i.e. deviations from the submanifold of coherent 
states) via an exact solution of the dynamics, which also 
appears novel to the literature. Finally we have applied 
our general results to coherent intertwiners in spin net¬ 
works as investigated in LQG. For expectation values of 
rotationally invariant operators (and these are the only 
ones contributing) one finds here a subleading correc¬ 
tion to the classical limit given in terms of universal (i.e. 
operator-independent) expansion coefficients which con¬ 
tain only geometric information about the network node. 


Let us now extend our results on the coherent-state 
expectation values of operator products to systems with 
N > 1 degrees of freedom. We start by two oscillatory 
degrees of freedom q a , p a with frequencies ui a , a € {1,2}. 
Iterating the arguments leading to Eq. (1361) one finds 


{a\AB\a) 


OO i-m —l— n 


m,n=0 


[Q+, 


a) 

(Al) 


with |a) = |aq) G> |a2) and 


Qa 


if, - -Pa 


(A2) 


Since Q i, Q 2 commute, the corresponding left arguments 
in the above nested commutators can be freely inter¬ 
changed such that 


blibl {c |[Q a , IQ,, A ]„] m | a) (a | [Qf, [Qf,B]J. 

= £ { a I [<9p m „(l)> \QPm. n ( 2))' • • [Qp mn (rn+n),A] ■ ■ ■]] | a) 




^P„(2)’ ' 


QPmn (m + n ) ’ ^ 


where the sum goes over all functions P mn '■ {1 ,,m+ value 1. Thus we arrive at 
n} —> {1,2} taking m times the value 2 and n times the 

I 


(A3) 


( a\AB\a) 


£tt£H [Qp*( 1)’ [QPn( 2)>- • ■ [Qp n (n),A\ •••]]!«) 

n =0 ' P n 



Q 


+ 

Pn(2)> 


Q~P n (n)’B 



(A4) 


where the second sum extends now over all functions 
P n : {1,..., n} —> {1,2}. Moreover, it is straightforward 
to see that the above expression also holds for an arbi¬ 
trary number N of oscillatory degrees of freedom with 


| a) = |cki) <E) ■ ■ • ® |ajv) and functions P n : {1,.. . ,n} —> 
{1,..., N}. In particular, the variance of an hermitian 
operator A can be expressed as 
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fcn 

(A Af =^ 2—^2 K“ \[QPn(l)’ [QPn{ 2)>‘" [Qp„(n),^] •■■]]!«) 

n—1 ' 


(A5) 


and the leading-order results for energy fluctuations are The counterpart of Eq. ED for two spins S \, S 2 reads 
given in Eqs. ED)-®. 


J 


2 S 2 2Si 


AAB\z) = Y, £ 


m —0 n =0 


(2 S 2 - m)\ (2Si - n)! 

m!(2S 2 )! n!(2Si)! 


-S~ 
h 2 ’ 


ft 5 ' 1 


1 0 + 
2 ’ 


_$+ n 

ft* 1 ’* 


(A6) 


with 1 2 :) = | 2 q) <g) l^). Due to the more complicated 
prefactors, a similarly compact form as in Eq. (IA4I) for 
the full expansion seems to be unachievable for spin sys¬ 
tems. The leading terms of energy fluctuations given in 
Eqs. mm, however, are again rather simple and al¬ 
low for an intuitive interpretation. 

Finally, combining both types of systems, the leading- 
order contribution to the fluctuation of an Hamiltonian 
depending on N oscillatory degrees of freedom and M 
spins reads 


N 


(A nf = 

a— 1 
M 

+£ 


Wa (dt.qa) 


{dtPa) 


O (h 2 ) 


b—1 L 


{hSb)2 ( 2 k {dtgb)2 + ° (J* 


(A7) 


Appendix B: The Normalization of Coherent 
Intertwiners and Related Integrals 


with 


(e z x s) x s 
\e z x si 


e z x s 
\e z x s| 


(B4) 


such that u, v, s, form an orthonormal system. If one had 
used the original operator m for generating coherent 
states the form of the vectors u, v would be less trans¬ 
parent. Now the normalization integral can be written 
as 


<*l P\*) 



(B5) 


with 


n n = nu a 


n% = nv a 


n n = ns a 


(B6) 


where n is the rotation axis occurring in Eqs. (EH),®. 
Taking into account the explicit form of the rotation ma¬ 
trix element 57 


S 


i'ipnS 


s) = 


ip . , . ip 

cos-1- in sin — 

2 2 


2 s 


(B7) 


In order to evaluate the normalization integral of co¬ 
herent intertwiners in the semiclassical regime, we shall 
use a slightly different version of SU(2) coherent states 
generated by 


V{Q,<p) = e“' pS V* s * (Bl) 

such that compared to Eq. mi one has (dropping again 
factors of h) 


V(d^)\S)=e^ s U(d,p)\S), (B2) 

i.e. the coherent states generated by the operators (USD 
and HD just differ by a phase factor which drops out 
from all expectation values. The operator (IB 1 1) fulfills 

V+SV = uS x + vS v + sS z (B3) 


elementary manipulations lead to 

■1 /»7T /»27T /»7t/2 

($|P|$) = -—77 / dd sind / dtp d^ sin 2 -0 

2^ J o Jo Jo 

N 

• ^ (^/cos^ + ins a sin 2Sa , (B8) 

? 7 =± a= 1 


where the cosine of ip is nonnegative in the entire integra¬ 
tion interval. Following RefAi we introduce p := n sin ip 
fulfilling 

d 3 p = sin 1 } sin 2 ip cos tpd'&dipdip (B9) 


such that 


($|P|$) 


J_ V /' —£^= e s "® 

27r2 tfr±Jp< 1 y/1 -p 2 


(BIO) 
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where 


S v (p) = ^ 2S a In (y \/1 - p 2 + zps a ) . (Bll) 

a—1 

In this form the integral can be evaluated via saddle point 
approximation to S v (p). As discussed in detail in Ref.—, 
provided that the classical closure relation (l95l) holds, the 
maximum of S v (p) occurs at p = 0 with 

N 

S+(0)=0 , S.(0) = 2niJ2Sa (B12) 

a— 1 

and since the latter sum must be integer for a nontriv¬ 
ial singlet space we have exp(S±(0))=l. The Hessian is 
given by (cf. Eq. (11011) ’) 


(B13) 

Extending now the integration domain in Eq. (IB 101) to 
the infinite space (as the integrand falls off rapidly), we 
are left with simple gaussian integrals leading to the re¬ 
sult (11001) where the leading first term was already ob¬ 
tained in Ref— while the subleading correction stems 
from expanding the square root in Eq. (IBlOj) . 

To compute the numerator of the coefficients HMD we 
consider 





■n<* 


e iil>n b Sb 



(B14) 


with V a = V(d a , fa)- With the help of the rotation ma¬ 
trix element¬ 


s'- 1 


0 iipnS 


S') = V2S ^cos + in z sin ^ 


25-1 


■(n x + in v ) sin ^ (B15) 


one derives 
S, 


iS x , e i ^” aSa 


iS y a , e ^ 3 ‘ 




S a ) = ~i ^cos ^ + ini sin 7 

•2 S a n y sin|-, 

o \ -I 'P , ■ z ■ 'P 

Sa) = 1 | cos — + m a sin — 


•2 S a n x a sin ^ , 


Sa) = 0. 


25-1 

(B16) 

25—1 

(B17) 

(B18) 


Now proceeding as before the two nontrivial expectation 
values can be formulated as 


(<t>\\iV a S x a V+,P}\$) 


-iSa 


E 

» 7 =±' 


d 3 p 


pV a 


o S r,{P) 


\Jl— p 2 rj^/l — p 2 + i 


'ipSa 


(B19) 


(<f>\\iv a sy a v a + ,p]\$) 

d 3 p 


_ iSg ST' 

7 r 2 

v=±' 


pUa 


1/1 - P 2 vV 1 - P 2 + ips a 


(B20) 


Performing again a saddle point approximation to the 
exponential and expanding the remaining integrand in 
quadratic order around p = 0 leads to 


(<S>\\iV a S x a V a + ,P]\<f>) 

(<!>\[iv a sy a v a + ,p]\<s>) 


, 1 ^n) 

\J 7 r det H 

{ulH^Sg) 
V 7 r det H 


(B21) 


(B22) 


and using Eq. (IB3I) along with elementary geometric re¬ 
lations it follows for the coefficients (l99l) 


3»(*) 


-S a (u a (v^H 1 S a ) - V a (UgH 1 S a )) 
SaSg X ^Vg {y a H S a ) 

+u a (u T a H~ 1 s a )y (B23) 


Finally, observing that u a , v a span the plane perpendic¬ 
ular to s a we obtain the result (I104D , and the numerator 
of Eq. (l99l) is given by Eq. (11031) . 
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